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ABSTRACT 



Aims. We study the statistical properties of the gravitational field generated by galaxy distribution 
observed by the Sloan Digital Sky Survey (DR7). We characterize the probability density function 
of gravitational force fluctuations and relate its limiting behaviors to the correlation properties of 
the underlying density field. In addition, we study whether the PDF converges to an asymptotic 
shape within sample volumes. 

Methods. We consider several volume-limited samples of the Sloan Digital Sky Survey and we 
compute the gravitational force probability density function (PDF). The gravitational force is 
computed in spheres of varying radius as is its PDF. 

Results. We find that (i) the PDF of the force displays features that can be understood in terms 
of galaxy two-point correlations and (ii) density fluctuations on the largest scales probed, i.e. 
r x 100 Mpc/h, still contribute significantly to the amplitude of the gravitational force. 
Conclusions. Our main conclusion is that fluctuations in the gravitational force field generated by 
galaxy structures are also relevant on scales ~ 100 Mpc/h. By assuming that the gravitational fluc- 
tuations in the galaxy distribution reflect those in the whole matter distribution, and that peculiar 
velocities and accelerations are simply correlated, we may conclude that large-scale fluctuations 
in the galaxy density field may be the source of the large-scale flows recently observed. 

Key words. Cosmology: observations; large-scale structure of Universe; 
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1. Introduction 

One of the main problems in studying the dynamics of a system of point masses concerns 
the analysis of the force acting on a test particle. This analysis can be achieved by determin- 
ing of the statistical properties of the (Newtonian) gravitational field generated by the given 
point-mass distribution. From this study in general one can derive several useful pieces of 
infor mation about the sta t istical and dynamical properties of the given point-mass distribu- 
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of 



tion (Chandrasekhar, 



Gabrielli , 



1943; 



2005 



Antonuccio-Deloau and Atrio-Barandela , 



Gabrielli et al. 
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, 



1992; 



Gabrielli et al., 



1999, 



20051) . It has been known since the pioneering article 



Chandrasekhar (1943) that, when particles are distributed without correlations (i.e., a Poisson 



distribution), the second moment of the gravitational force diverges. Certain spatial correlations be 



tween particles may cause the divergence, for instance, of the first moment (iGabrielli et al 



1999). 
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For this reason, instead of computing the moments of the gravitational force, it is more appropriate 
to determine its probability density function (PDF). In particular, to understand the relevant features 
of the gravitational force fluctuations, it is interesting to study the limiting behaviors of the PDF in 
both strong and weak fields. It can be shown (see below) that the former is closely related to the 
small scale two-point correlation properties of the distribution. On the other hand, the weak field 
behavior is determined, in a non-trivial way, by how rapidly the distribution becomes isotropic on 
large scales. There fore, this latter behavior is connec ted to the higher order correlation properties 



of the distribution dGabrielli et all 



1999, 



2005, 



2006) 



Our first aim, in this paper, is to consider the determination of the PDF of the gravitational 
force generated by the observed spatial distribution of galaxies. Given that the galaxy distribution 
is described as a stochastic point process, the statistical properties of the gravitational field on an 
arbitrary point is itself a stochastic quantity. The properties of these two stochastic processes are 
closely related. In particular, the force PDF can be simply regarded as a particular statistical char- 
acterization of the galaxy distribution containing some useful statistical information about higher 
order galaxy correlations. 

To connect the force PDF of the galaxies with that of the whole mass distribution, we have 
to make an hypothesis about the relation between luminous and dark matter. The most simple 
assumption is that the galaxies trace the peaks of the underlying matter density field, i.e., that the 
galaxy density field is linearly biased with respect to the dark matter density field. If this is so, 
then the fluctuations in the gravitational force in the galaxy distribution also reflect those in the full 
matter distribution. 

This is actually the standard hypothesis that is usually assumed in studies that compare of the 
dipole induced by the gravitational influence of structures in our local Universe with the observed 
CMB dipole [j. Nowadays, several groups have sho wn that structures responsible for the dipole 



must be located on scales larger than ~ 150 Mpc/h (Mail er et al 



2003). However, because of the 



sparseness of data at very large distances and the incompleteness of galaxy surveys, ther e is still no 



2008). 



consensus about the depth of the convergence for the CMB dipole (ILavaux et all \ 

The actual dipole caused by structures around us corresponds to a determination of the force 
around a single point of the galaxy distribution (our galaxy). On the other hand, the PDF is de- 
termined by computing the gravitational force on each galaxy contained in a given catalog, and 
by repeating this measurement for all measured galaxies. With respect to the single determination, 
the PDF has several advantages: (i) it considers complete samples extracted from galaxy redshift 
surveys (i.e., volume-limited samples) without making any additional assumption about the incom- 
pleteness of the catalogs, (ii) We do not have to consider the sky regions where observations were 
not done and use weighting schemes to correct for this incompleteness. We indeed compute the 
force only in spherical volumes fully contained in our samples, (iii) By considering each galaxy as 
a center, we determine the force generated by all galaxies contained in a spherical volume around 
it. We can therefore evaluate the statistical significance that the determination of the force from a 
single point (i.e., the dipole) is (or is not) found to converge to an asymptotic value within a certain 
distance. 



1 This shows that our Local Group moves with a velocity of 627 ± 22 km/ s towards the direction / = 276° 
and b = 306° in Galactic coordinates. 
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It is then interesting to consider the relation between the gravitational force field g and the 
peculiar velocity field v. In the linear regime of low-amplitude density fluctuations, it is well- 
known that the velocity is proportional to the gravitational acceleration (IPeeblesL Il980l) . As our 



observations are limited to a range of scales in which galaxy fluctuations are still large, we can 
only speculate about the relation between the gravitational force and the peculiar velocity fields 
by assuming that there is a direct proportionality even in the non-linear regime. In this respect we 
note that when fluctuations are large and t he universe is curva ture-dominated (i.e., with negligible 
matter content) one still finds that v oc g ([Joyce et al. 1 120001) . The relation between g and v has 
acquired interest because of the increasing number of observations of peculiar velocities on large 
scales. 

Peculiar velocities provide important dynamical information because they are related to the 
large-scale matter distribution. By studying their local amplitudes and directions, these velocities 
allow us, in principle, to probe deeper, or hidden, parts of the Universe. The peculiar velocities are 
indeed directly sensitive to the total matter content, inferred from its gravitational effects, and not 
only to the luminous matter distribution. However, their direct observation by means if distance 
measurements remains a difficult task. There have been observations of larg e-scale galaxy coher- 

J2OO9I) 



Watkins et al. 



ent motions that are at odds with standard cosmological models. For instance 

measured in a Gaussian window of diameter lOOMpc/h, a coherent bulk motion of 407 + 81 km/s 

which conflicts with the prediction of « 200 km/s at the 2<x le vel of LCDM. Other independent re 



suits confirm the presence of unexpectedly l arge b ulk motions dLavaux et al 



2008 



2009). In particular, 



Kaslinski et al 



2008; 



Kasli nski et all 



(2008) measured peculiar velocities of galaxy clusters 
by studying, in a large statistical sample of clusters, the fluctuations in the cosmic microwave back- 
ground (CMB) generated by the kinematic Sunyaev-Zeldovich (KSZ) effect, i.e., caused by the 
scattering of microwave photons by the hot X-ray emitting gas inside clusters. They measured a 
bulk flow of 600-1000 km/s on scales > 300 Mpc/h (the limit of the considered clusters catalog) 
which is more than 10 times larger than expected by the LCDM model. In addition, they concluded 
that the coherence len gth of the measured bulk flow shows no signs of convergence out to scales. 
(2009), by applying the KSZ method to a larger cluster sample, showed that the 



Kaslinski et al 



flow is consistent with approxima tely constant v elocity out to at least ^ 800 Mpc/h. Although 
these results remain controversial JKeislen , |2009), measurements of peculiar velocities have very 



interesting potentials a nd future surveys are expected to prov ide more precise data of bulk flow 



motion on large scales dSmith et al 



2004; 



Frieman et all 



2008). Our aim here is to verify whether, 



based on the assumption that v oc g 7 fluctuations in the force distribution should be compatible 
with the large-scale bulk flows observed. If this were so, this would indicate that local galaxy struc- 
tures may represent an important contribution to the peculiar velocities that could avoid us having 
to invoke exotic physi cs, as the effect of far- away pre-inflationary inhomogeneities, to explain the 



large-scale bulk flows (iKaslinski et al 



2008). 

In this paper, we focus on the following two questions: (i) What is the shape of the PDF and 
how is this connected to the underlying fluctuations in the galaxy density field ? (ii) Is there a 
convergence of the PDF to an asymptotic shape that is independent of the size of the volume (i.e., 
the sphere radius r) in which this is computed. Based on the assumption that light traces mass, we 
are unable to determine the direction and amplitude of the dipole around us, thus verify whether 
this is the same as the CMB dipole. However, we can characterize the properties of gravitational 
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fluctuations generated by galaxy structures in a complete statistical way (up to the maximum scale 
R s allowed by the geometry of available samples); in particular, we can determine whether in the 
volumes considered the gravitational force has reached, from a statistical point of view, saturation 
on scales rj < R s . In this case, fluctuations on scales larger than rf do not contribute anymore, 
statistically, to the amplitude of the gravitational force. That is, large-scale contributions to the 
force from different directions tend to compensate each other. The scale rj would thus be the scale 
beyond which large scale flows would start to decay in amplitude. On the other hand, if we were 
able to place only a lower limit on the scale rf, i.e. if ry > R s without a convergence of PDF to a 
stable shape within the sample volume, this would corroborate the conjecture that large-scale flows 
can be sourced by fluctuations in the density field that statistically break spherical symmetry. In this 
case, larger and larger volumes would continue to contribute to the gravitational force on an average 
galaxy. Therefore, if the measurements of large-scale bulk flows are unaffected by major statistical 
and systematical effects, but are generated instead by galaxy density inhomogeneities, we would 
expect not to see a convergence of the PDF to a stable shape inside the considered samples because 
we are limited to a volume of radius R s w 100 Mpc/h. As discussed above, while the information 
on the force PDF of galaxies can be derived in straightforward way from the data, its relation to the 
peculiar velocity field is based on a series of assumptions. Nevertheless, it is interesting to consider 
whether when galaxies trace mass and peculiar velocities are parallel to accelerations, it is possible 
to learn something about the relation between the force and peculiar velocity fields. 

The paper is organized as follows. In Sectf2] we briefly review some basic properties of the 
gravitational force field in stochastic point distributions. We then describe the selection of the 
galaxy samples in Sectj3]and the methods (Sectj4]i used to estimate the relevant statistical quanti- 
ties in the data. In Sect[5] we discuss our main results and in Sect|6]we consider the PDF of the 
gravitational force in mock galaxy catalogs generated from cosmological N-body simulations. We 
finally draw our main conclusions in Sect|7] 

2. Gravitational force from stochastic point distributions 

The value of the gravitational force acting on a particle belonging to a given point distribution is 
determined by both the immediate neighborhood of the particle itself and the large-scale proper- 
ties of the system. This situation reflects two peculiar features of the gravitational force, that it 
is divergent at small separations and is slowly, i.e., in accordance with a power-law, decaying at 
large distances. The stochastic nature of a given point distribution generates a certain PDF P(F) 
of the gravitational force F. In general, one would like to calculate P(F) from the statistical prop- 
erties of the underlying particle distribution. There are a few examples of particle distributions 



with analytically calculable P(F). IChandrasekharl d 1943b was the first to compute the PDF of the 



gravitational force for an uncorrected distribution of points, i.e. the Poisson distribution. In this 
situation, the PDF of the force is found to be given exactly by the Holtzmark distribution (HD), 
which is a three-dimensional fat (or heavy) tailed Levy distribution. Approximated generalizations 
of this appro ach can be found for more complex particle systems obtained by pertu rbing a Poisson 



2006) 



2004) or perturbing a perfect cubic lattice of particles (Gabrielli et al 



distribution (Gabriell i et al. 

Vald (1994) developed a functional integral approach for evaluating the stochastic proper- 



ties of vectorial additive random fields generated by a variable number of point sources obeying 



Sylos Labini: Gravitational fluctuations of the galaxy distribution 



inhomogeneous Poisson statistics. In this way he derived the generalized Holtzmark distribution 
(GHD), which we discuss in more detail below. 

The gravitational force for unit mass due to the N particles contained in a sphere C(r, jc,-) of 
radius r centered on the i' h particle of a given distribution is 



Fi(r) 



Z 



Gm 



7 !*,■-*/ 



(l) 



where G is the gravitational constant, nij is the f h particle mass, and N is the number of particles 
in a given volume V. In what follows, we take G - 1 and we consider particles with the same mass, 
so that we assume that m,- =1 V i. 

If the particle distribution is statistically stationary, i.e. it is invariant under space rotations and 
translations, the direction of F has equal probability in each direction, i.e. P(F) depends only on 
the force modulus F - \F\. For this reason, we can consider the PDF W(F) of the force modulus 
F, which is simply related to the three dimensional PDF by W(F) - 4jtF 2 P(F) . 

The statistical properti es of a density field can be characterized by measuring the average con- 



ditional density defined as dSvlos Labini et al. , 



2009a) 



(n(r)) p 



1 N 

/=1 



(2) 



where rtj(r) is the density measured in a sphere of radius r centered on the i th particle, and the 
second equality holds only for the case in which there are power-law correlations in a certain range 
of scales. In Eqf2] the constant A is related to the average nearest neighbor (NN) distance A « A~ l ^ D 
and y — 3 - D is the correlation exponent^ . When there are power-law correlations on small scales 
one can readily derive an exp ression for the force PDF due only to a NN particle by using the 



identity (IGabrielli et al. , 



2005) 



W m (F)\dF\ = u(r)\dr\ 

where 

io(f) = \nAr D - 1 



4jtA , 

cx -P\~— r 



(3) 



(4) 



is the NN PDF for a particle system obeying to Eq|2"l 3 , For D = 3, Eq|4]reduces to the well-known 

1943). Thus, from Eqf3]|4]we can derive the 



expression for the Poisson NN PDF dChandrasekhar 
NN force PDF 



D 



W nn (F;D)=jF A 



' F - exp 



-m" 



(5) 



\2/D 



where F\ - (^f A) . For D — 3, Eqf5]reduces to the Poisson case dChandrasekhar 
asymptotic behavior in the strong field limit is 



1943). The 



D 



\xmW nn {F;D)^-F u ^F-- 



D/2j 



(6) 



F^m 2 

It is simple to show that the average force diverges for D > 2 and that the average square force 
diverges VZ) < 3. For this reason, as mentioned in the introduction, we do not directly consider the 
study of these moments, while we focus on the shape of the PDF. 



Note that D is no t the space dimension, but it can be interpreted as the fractal dimension of the object 



i Gabrielli et al. . 



2005). 



More precisely, Eq|4]is found when higher order correlations are neglected. 
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Valdl (119941) derived a simple expression for the GHD, which holds when there are only 



p ower-law two-point co rrelations as in Eqf2] and no higher order ones (see also the discussion 



Gabrielli et al. 



( 1999)). The GHD can be written as 



W(F; D) 



H(J3-D) 



F 
where /3 = F/Fq, 



2 r°° 

H(B; D) = — exp (-(x/B) Dn ) x %m{x)dx , 
np Jo v ; 



(7) 



(8) 



the constant F^ is given by 
( 
(D + 2)r £ (f)sin(f)J 



2/D 



(9) 



and the symbol Eg represents the Euler Gamma function. For D — 3, Eqs|7]|9]reduces to the HD. 
It is easy to show that in the strong field limit W(F; D) — > W„ n {F; D), so that the fat tail is always 
present and due to NN correlations. 

The PDF of the vector F, i.e. P(F), must be symmetrical about zero, because the distribution 
is isotropic and the average of the vector force must be zero. Thus P(0) = const., so that we have 



lim W(F) oc F 2 P(0) oc F 2 



(10) 



which is indeed the case for both the HD and the GHD. Hence the limiting behaviors of the GHD 
(and of the HD for D = 3) are 



W(F) oc F- (D+2)/2 for F 
W(F) oc F 2 for F -> . 



(11) 



The strong field limit is determined by the small-scale correlation properties, while the weak field 
limit is very general and related to the stationarity of the distribution. We focus on these tails in the 
analysis of the galaxy samples. 

We emphasize that Eqf2] may describe an highly inhomo geneous structure c haracterized by 



2005). In this a situ- 



large fluctuations and non-trivial higher order correlations (Gabrielli e t all 

ation neglecting higher order correlations is a very crude approximation. However, only in this 

limit is one able to calculate the force PDF. By taking into account higher order correlations, 



one may compute only few moments of the force (IGabrielli et al 



1999, 



2005). This means that, 



while the limiting behaviors given by EqQT]are expected to be satisfied as long as Eqf2] holds and 
the particle distribution is stationary, the exact location of the peak, marking the transition from 
W(F) oc /M£>+2)/2 to y/(j?^ x p2^ j s poorly constrained from a theoretical point of view. 

We note that Eqf2]may also describe a structure that has varying correlations as a function of 
scale. In this case, D = D(r) < 3 where D — 3 when the distribution becomes spatially uniform and 
the average conditional density coincides with the uncondition al density, modulo s mall-amplitude 
correlations described by the reduced correlation function £(r) (Gab rielli et alll2005l) . Only in this 
situation may one describe the properti es of the force PDF in terms of either the f (r) fun ction or its 



Fourier transform, the power spectrum (IGabrielli et all 



2004, 



2006, 



2005; 



Joyce, 



2007). 
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VL sample 


"min 


"max 


w*-min 


'-''-max 


N 


A 


VL1 


70 


450 


-20.8 


-21.8 


112860 


1.5 


VL2 


50 


200 


-18.9 


-21.1 


73810 


0.6 


VL3 


200 


600 


-21.6 


-22.8 


51697 


2.6 



Table 1. Main properties of the obtained VL samples with K-corrections and without E- 
corrections: R mm , R max (in Mpc/h) are the chosen limits for the metric distance; M m ,„, M max define 
the interval for the absolute magnitude in each sample. Finally, N is the number of galaxies in the 
sample and A is the nearest neighbor distance (in Mpc/h). 



3. The data 



We extracted a volum e-limited (VL) s ample constructed from t he data rel e ase 7 (DR7) 



(LAbazaiianet al. , 



2009) of the SDSS (see 



Svlos Labini et al. 



(2009a); 



Antal et al. 



(2009) for de- 



tails). To compute the metric distance, we used the standard cosmological parameters Q^ = 0.3 
and Qa = 0.7. We considered a contiguous sky area with almost uniform coverage and we applied 
K-c orrections to calculate abs olute magnitude M, in the r band (corrected for Galactic absorption) 
(2009a). 



as in 



Svlos Labini et al 



We selected three samples (see Table [TJ. The sample VL1, containing about one fifth of all 
galaxies in DR7, has relatively large spatial extensions and small spread in galaxy luminosity. On 
the other hand, the sample VL2 contains fainter galaxies and covers a smaller volume, while VL3 
contains brighter galaxies and is more extended in radial distance. 

In Table Q] we also report the average distance between nearest neighboring galaxies: this pro- 
vides a lower cut-off for the computation of real-space statistical properties. It is also directly 
related to the l arge-field tail of the P DF of the force. The fainter the galaxies, clearly, the smaller 



the value of A (Gab rielli et al. , 



2005). 

We do not apply any correction forredshift distortion,thus our results are given in redshift space. 
We expect the peculiar velocities to affect the small-scale properties, where their redshift distortions 
are relatively important, of the statistical quantities we measured but to leave conclusions on large 
scales unchanged. 



4. Methods 

The average conditional density in Eqf2] is computed by evaluating, on the scale r, the average 
of the N(r) density determinations over the sample points whose minimal distance from sample 
boundaries is smaller than r. The force F (Eq[TJi is computed in a sphere of radius r centered on 
a galaxy with the same constraint as above, so that only comple te spheres are considered . This 



2009a). We 



implies that the number of determinations depends on the scale r dSvlos Labini et all 
limit the analysis to N(r) > 10 4 . We note that the force PDF we compute here is intrinsically 
defined for a discrete and stochastic point distribution. As long as these points (galaxies) display 
non-trivial spatial correlations, the force PDF deviates from the Holtzmark distribution. 
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5. Results 

We find that the average conditional density for the VL1 sample (see FigQ] inset panel) behaves 
for r < 20 Mpc/h as 



(n(r)) p = 0.036 x r 



-0.9 



while for 20 < r < 80 Mpc/h we get {n(r)) p = 0.0058 x r^ lb (lAntal et al. . 



200' 



(12) 



This be- 



havior implies that a crossover toward homogeneity not being found up to ~ 80 Mpc/h. This im- 
plies that fluctuations are still large when filtered on scales ~ 80 Mpc/h and that the standard 
two-point correlation function ^(r) results, which are biased by major systematic effects in these 
sa mples, may undere s timate the amplitude of fluctuations on scales > 10 Mpc/h (see the discussion 



Svlos Labini et al. 



(I2009aldln . 

In other two samples considered, VL2 and VL3, we found similar behaviors (see Figsj2j3}: (i) 
the PDF of the force does not converge to an asymptotic shape up to ~ 100 Mpc/h and (ii) the 
conditional density correspondingly does not exhibit a well defined trend toward homogenization. 




Fig. 1. PDF of the force for the VL1 computed one different scales r. The solid line is EqQwith 
the parameters from Eqfl~2] In the inset panel, the behavior of the average conditional density is 
plotted in addition to its best fits for r < 20 Mpc/h and for 20 < r < 80 Mpc/h. 



In FigQ] the PDF of the force for the SDSS-VL1 is shown at different scales r. We note that the 
large field tail is stable for all sphere radii considered. It is characterized by a power-law behavior 
of the type W(F) ~ F~ a with a - -2.1. From Eq|5l we infer that D - 2.2, thus in agreement with 
the value found for the conditional density on small scales (i.e., EqlTZl). Thus the validity of the 
approximation at large F values can be simply explained by very small-scale considerations. For 
r — 30 Mpc/h, the GHD obtained from Eqs|7]|9]with the values of A,D measured from the best 



4 The relation between the conditional density and the standard two-point correlation f unction, and in 
particular the relation between their respective exponents in a finite sample is discussed in 



Gabrielli et al 



(2005) (see pg.240-242). 
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Fig. 2. The same of FigQ]but for the sample VL2 




Fig. 3. The same of FigQ]but for the sample VL3 



fit of the conditional density for r < 20 Mpc/h, fits the data extremely well. (i.e. Eq[T2l. In this 
case, there are no free parameters to be adjusted. In addition, we note that the weak force tail is 
proportional to F 2 , in agreement with EqfTTl 

For r - 90 Mpc/h, the same fit with the GHD, for small F values, is not as good as for the pre- 
vious case, while it behaves again as F 2 in agreement with EqQj] This is because the contribution 
to the gravitational field from the largest scales in this sample are non-negligible, and in addition 
the average conditional density displays a different slope for r > 20 Mpc/h than for r < 20 Mpc/h. 
This implies that the GHD distribution is only an approximation to the real PDF. It is indeed clear 
that the simplifying assumptions used to derive Eq|7]|9j and in particular the hypothesis that there 
are no n-point correlations, are too naive with respect to the complexity of the true galaxy struc- 



10 
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ture. The value of the field at which the turnover of the PDF occurs, depends on the sphere radius r 
determining the cut-off to the force calculation in the data. This implies that up to the largest sphere 
radius considered we do not find a convergence of the force PDF to an asymptotic behavior. The 
PDF peak instead shifts toward stronger field values at larger r, implying that greater and greater 
contributions are given by the large scale fluctuations. We note that there is a factor of three be- 
tween the PDF peak location at r — 10 Mpc/h and at r — 90 Mpc/h while the peak of the PDF 
increases by 20% from r — 50 Mpc/h to r — 90 Mpc/h. It is surely smaller than from r — 10 Mpc/h 
to r — 50 Mpc/h but this implies that a non-negligible contribution to the force field is caused by 
large-scale fluctuations. 



6. Comparison with mock galaxy catalogs 

To compare the results obtained in the data with the theoretical predictions, we considered 
a semi-analytic gala xy catalog constructed from the Millennium LCDM N-body simulation 
dSpringel et all 120051) . To construct mock samples corresponding to SDSS VL samples, we used 
the full version of the cat alog in the usriz filte r system. The catalog contains about 9 million galax- 



*>B 



ies in a 500 Mpc/h cube (ICroton et all 2006) |j. We use the absolute magnitudes in r filter used in 



the SDSS case, to construct the mock samples with the same limits in absolute magnitude as for 
the SDSS VL samples with K-corrections. As for the real data, we considered the redshift-space 



Svlos Labini et al. 



distribution (more details on the construction of mock samples can be found in 
d2009al) ). The average distance between the nearest points in the mock galaxy sample for which we 
present results is A ~ 1 .6Mpc/h, i.e. about that of the sample VL1 of SDSS. We note however that 
in other mock samples we found very similar results to the ones discussed here. 

The main features of the considered statistics in the mock galaxy samples are (see PigSJ (i) 
beyond a few tens of Mpc the PDF stabilizes to an asymptotic (sample-independent) shape and 
(ii) the conditional density exhibit a well-defined flattening on large enough scales. Both these 
behaviors differ from those found in the real data. 

To conclude this section, we note that in LCDM models, on large enough scales, one expects 
the force PDF to be sample-size independent. In general, if the variance of the force is finite, the 
force PDF is necessarily w ell defined. F rom the definition of the force, and after transforming into 



Fourier space we find that dJovce , 



2007) 



(F 2 )oc jd 3 kk~ 2 P(k), 

where P(k) is the power spectrum of the density field. Therefore, a sufficient conditio n for the forc e 
to be defined in the usual thermodynamic limit is just that k~ 2 P(k) be integrable dJovce , 



2007) 



Since we are interested here in only the possible divergences in this quantity caused by the long 
distance behavior of the fluctuations, it is therefore sufficient that we assume that lim*_>o kP(k) = 0. 
Thus it follows that the force is well defined in the infinite volume limit if P(k) diverges more 
slowly at small k than k~ l . This condition is satisfied by st andard cosmological density field s for 



which P(k) ~ k on large enough scales, i.e. r > 100 Mpc/h dSvlos Labini and Vasilvev , 



2008). 



5 See http://www.mpa-garching.mpg.de/galform/agnpaper/ for semi-analytic galaxy data files 
and description, and see http://www.mpa-garching.mpg.de/millennium/ for information on 
Millennium LCDM N-body simulation. 



Sylos Labini: Gravitational fluctuations of the galaxy distribution 



11 





1 


'L>,-- 


1 i.ii 





- 




I ft 


A"^! 








: 


10° 




— r=20 
•—r=60 
— . r=120 








to' 1 


rjN Ip Jr 








- 




_ N 1 i 


"' 


r 

10° 


io 1 io ! 






" 








1 ' ' 


III 1 1 












to" 2 


io" 1 
- A 

v 10 - : 


- 




const. 




Is 




" 








3 




1 . . 


ni N.t,"n 




,uL 


- 






. ( 




i 


, , 



Fig. 4. PDF of the force for the mock sample computed on different scales r. In the inset panel, we 
plot the behavior of the average conditional density. 



7. Conclusions 

The gravitational force probability density (PDF) in finite volumes represents an interesting char- 
acterization of the statistical properties of a given point distribution which is directly related to the 
statistical properties of density fluctuations and their dynamical effects on large scales. By analyz- 
ing only the statistical properties of the force PDF we have found that this PDF is in general related 
to two-point and higher order correlations. We have measured the PDF in several galaxy samples of 
the SDSS finding that (i) its limiting behaviors can be simply inferred from the galaxy correlation 
properties and that (ii) the PDF does not achieve an asymptotic, sample-size-independent, shape 
within the volumes considered, i.e. ~ 100 Mpc/h. These results agree with the studies of the local 
dipole due to structures, and in particular those that do not find a convergence up to ~ 100 Mpc/h. 

The relation between the force PDF and the dynamics, i.e. between gravitational accelera- 
tion and peculiar velocity fields requires several comments. First, we have to assume that galaxies 
trace the underlying matter density field, at least on large enough scales (i.e., that biasing is lin- 
ear). This is a reasonable assumption usually employed in the studies of the dipole due to struc- 



tures ( Mailer et al. 



2003; 



Lavaux et al 



20081) . It is then necessary to assume that peculiar veloc- 



ities are paralle l to gravitational accelerations. This can be verified by linear perturbation theory 



(Peebles, 



dJovce et al. 



1980) or when fluctuations are large but the universe is in a curvature-dominated phase 



2000). By using these assumptions, we conclude that the observed values of the bulk 



flows are compatible with the gravitational force fluctuation field from galaxy structures. 

In other words, these results support the conjecture that large-scale structures may be the source 
of the observed large-scale flows, which are unexpected in the standard LCDM scenario. We note 
that even the correlation properties of the galaxy density field in the SDSS samples are also at odds 



with the standard LCDM predictions (ISvlos Labini et al. , 



2009 



HH 



e). 
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The results of this paper corroborate that the observed large-scale flows are due to large-scale 
fluctuations in the galaxy density field. A more quantitative conclusion must consider the dynamical 
history of structure formation and is beyond the scope of this paper. 
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